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Abstract 

Interconnection and damping assignment passivity-based control (IDA-PBC) is an 
excellent method to stabilize mechanical systems in the Hamiltonian formalism. In this 
paper, several improvements are made on the IDA-PBC method. The skew-symmetric 
interconnection submatrix in the conventional form of IDA-PBC is shown to have 
some redundancy for systems with the number of degrees of freedom greater than 
two, containing unnecessary components that do not contribute to the dynamics. To 
completely remove this redundancy, use of quadratic gyroscopic forces is proposed 
in place of the skew-symmetric interconnection submatrix. Reduction of the number 
of the matching partial differential equations in IDA-PBC, and simplification of the 
structure of the matching partial differential equations are achieved by eliminating the 
gyroscopic force from the matching partial differential equations. In addition, easily 
verifiable criteria arc provided for Lyapunov/exponential stabilizability by IDA-PBC 
for all linear controlled Hamiltonian systems with arbitrary degree of underactuation 
and for all nonlinear controlled Hamiltonian systems with one degree of underactuation. 
A general design procedure for IDA-PBC is given and illustrated with two examples. 
This paper renders the IDA-PBC method as powerful as the controlled Lagrangian 
method that is a Lagrangian counterpart of IDA-PBC. 



1 Introduction 



Mechanical systems are ubiquitous in nature and engineering, and there have been many 
studies on modeling motions of insects, animals, fish and humans in the framework of 
mechanics and control. For example, a model for clock-actuated legged locomotion of hu- 
man and biologically-inspired robots is proposed and analyzed in [20J. A control-theoretic 
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strategy for human walking gait assistance is suggested using the biped model to lessen the 
perceived weight of a patient's center of mass through a robotic angle-foot orthosis with 
one actuated degree-of- freedom |15j . A time-scaling control law is developed and applied 
to two passive-dynamic bipeds: the compass-gait biped and a simple biped with torso [16] . 
In all of these works, one of the main objectives is to study stability and stabilization of 
the motion of systems. 

The energy shaping method stands out among several methods for stabilization of 
mechanical systems since it preserves the mechanical structure, provides a systematic pro- 
cedure for constructing control laws, and yields a large region of stability for closed- loop 
systems. The idea of the energy shaping method is simple. Given an unstable mechan- 
ical system, one transforms it via feedback to a stable mechanical system whose total 
energy function obtains a minimum value at the equilibrium of interest. In this process of 
transformation, the original total energy function with a saddle-type critical point at the 
equilibrium point gets transformed to a new total energy function with a minimum value 
at the equilibrium point. In order to find such a new mechanical system with a stable 
energy function, one has to solve partial differential equations (PDEs) for the mass matrix 
and the potential function of the new mechanical system. The PDEs for the new mass 
matrix are called kinetic matching conditions or kinetic matching PDEs, and the PDEs for 
the new potential function are called potential matching conditions or potential matching 
PDEs. Hence, understanding the structure of the matching PDEs and their solvability is 
important. 

The idea of potential energy shaping for stabilization of a mechanical system dates 
back to [2J. The notion of kinetic energy shaping first appeared in [3]. The idea of 
total energy shaping was introduced in [U [6], and has then been actively developed 
El [TOl HI US HH [21]. There are two approaches to energy shaping: the La- 
grangian approach and the Hamiltonian approach. The energy shaping method is called 
the method of controlled Lagrangians on the Lagrangian side and the method of intercon- 
nection and damping assignment passivity-based control (IDA-PBC) on the Hamiltonian 
side. It is proven that the two approaches are equivalent [8]. Chang [9] has then improved 
the Lagrangian method by completely characterizing quadratic gyroscopic forces, reducing 
the number of matching conditions for energy shaping, and finding necessary and sufficient 
conditions for stabilizability by energy shaping for the class of all linear mechanical systems 
with arbitrary degree of underactuation and the class of all nonlinear mechanical systems 
with underactuation degree one. In contrast to these developments on the Lagrangian side, 
there has been a lag in the development on the Hamiltonian side. 

In this paper, we make improvement of the method of IDA-PBC in several ways so 
that the IDA-PBC method becomes as strong as its Lagrangian counterpart. First, we 
propose a new form of IDA-PBC by introducing a quadratic gyroscopic force. The skew- 
symmetric interconnection matrix in the conventional form of IDA-PBC [Hill], is shown 
to have unnecessary components that do not appear in the system dynamics with the num- 
ber of degrees of freedom greater than two. To remove the unnecessary components, use 
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of quadratic gyroscopic force is proposed in place of the skew-symmetric interconnection 
matrix. Second, the kinetic matching conditions get decomposed into two sets: one con- 
taining no components of the gyroscopic force and the other containing components of the 
gyroscopic force. The first set constitutes a new reduced set of kinetic matching PDEs 
and the second is used to algebraically define the gyroscopic force. This decomposition 
substantially reduces the number of kinetic matching PDEs. For example, if the degree of 
underactuation is one, there is only one kinetic matching PDE, irrespective of the num- 
ber of degrees of freedom. Moreover, when the co-distribution spanned by actuation (or 
control) co-vector fields is integrable, the new reduced matching PDEs contain a smaller 
number of entries of the unknown mass matrix, which was never discovered in the past. 
Third, we provide necessary and sufficient conditions for stabilizability by IDA-PBC for 
the class of all linear mechanical systems with an arbitrary degree of underactuation and 
for the class of all nonlinear mechanical systems with one degree of underactuation. These 
conditions can be easily verified in advance before solving the matching PDEs. Fourth, 
a step-by-step synthesis procedure with IDA-PBC is provided and the main results are 
illustrated with two examples. 

2 Main Results 

We review some basic notions on tensors, derive important lemmas, present a new form 
of IDA-PBC and compare it with the conventional form. Although both forms are shown 
to be theoretically equivalent for use in IDA-PBC, the new one has the advantage that 
it contains fewer components. We then exclusively work with the new form. We derive 
matching conditions for energy shaping and decompose them into two parts: one without 
the gyroscopic term and the other with gyroscopic term. This decomposition allows to 
reduce the number of matching PDEs. We discuss conditions for stabilizability by IDA- 
PBC that can be verified easily, and illustrate the results with two examples. 

2.1 Important Lemmas 

Let V be an n-dimensional real vector space and V* its dual space. Let {ei, . . . ,e n } be a 
basis of V, and {e 1 , . . . , e n } its dual basis such that (e i ,ej) = <5j, where (, } is the canonical 
pairing between dual vectors and vectors, and 5j is the Kronecker delta. The tensor product 
of two vector spaces V and W is denoted by V ® W, and each element of V ® W is a linear 
combination of elements of the form v <S> w, where v € V and w e W. The r-fold tensor 
product of a vector space V is denoted by V® r or ® r V. An (r, s)-tensor T on V is written 
as 

T = T^;^e il «•••«(;. Z( J Z---Z e js , 

where the Einstein summation convention is enforced. The contraction of an (r, 0)-tensor 
S = S^'^e^ <8> ••• <g> e ir and a (0, s)-tensor T = Tj v ..j s e^ <S> ■ ■ ■ <2> e^ 3 with r < s, is defined and 
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denoted by 

S j T = S i ^T ll ... irlr+v .. i /^ ® ... ® e *\ 
One can identify each (0, s)-tensor T with a multi- linear map T : V x ••• x V -»■ R defined by 

T(«i,— ,v s ) = «i ® — ®u s jT 

= v s j ■ ■ • j v\ j r. 

Define the symmetrization operator Sym : V® r -> y® r by 

Sym(«i ® -®v r ) = — £ ^ ct(1) ® ••• ® t)^), 

where S r is the symmetric group of {l,...,r} and it is understood that Sym is linearly 
extended to V® r . 

A (0, 3)-tensor C = Cy&e 1 ® e- 7 ® e k on V is said to be gyroscopic if it satisfies 

C(u,v,w) = C(v,u,w) \/u,v,weV, (1) 
C(u,v,w) + C(v,w,u) + C(w,u,v) = Vii,u,wey. (2) 

In coordinates, JJJ) and (J2J) are equivalent, respectively, to the following 

Cijk - Cjik, (3) 
Cijk + Cjki + Ckij = (4) 

for all i,j,k. Define 

C(V) = {<7 e V*® 3 | properties © and © hold} 



{Cij k e l ® e j ® e fc | properties © and © hold}. 



Define 



B(V) = {5 e F* 03 | B(u,v,w) = -B(u,w,v)Vu,v,w e y} 
= {Sijfce* ® e J ® e k \ B ijk = -B ikj }. 

Lemma 1. Let n = dimF. Then, dim£>(V) = n ^""^ and dimC(V) = n ( n ^ . 

Proof. Due to the skew-symmetry property B^k = -Bikj, we have dim£>(V) = !ll|zil To 
compute dimC(y), we consider three cases: 1) all three identical indices; 2) two identical 
indices and one distinct index; and 3) three distinct indices. For case 1, we only have 
Cm = by ©. For case 2, we have C^j = Cku = ~\Cak for all i + k by ([5]) and (JH). 
Hence, the number of independent components in this case is n{n - 1). For case 3, the 
number of independent components is n ( n-1 )( n 2 ) by p anc [ Q. Hence, d\mC{V) = 
+ n(n - 1) + "("- X K"- 2 ) a 2^=12. ' □ 
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Lemma 2. Let tp : B(V) -> V* 03 be a linear map defined as 

r(i '>;,/,■>' ® e j ® e k ) = -(B ijk + B jik )e i ® e> ® e k (5) 

for all Bij k e { ® e j ® e fe e B(T^). T/ien, V(#(^)) = C(V), i.e., i/ie ima#e of B{V) under the 
map if) equals C(V). Moreover, dimker^ = n ( n 1 )( n 2 ) ; where n = dimV. 

Proof. First, we show that tp{B{V)) c C(V). Choose any B = By^e' ® e- 7 ® e fc e 
Let C = ip(B). In coordinates, Cjjfc = \{Bij k + Bj ik ), where C = C^e 1 ® e 1 ® e fc . Then, 
Cyfc = \{B%ik + -Bjifc) = \{Bjik + Bij k ) = Cj ik , satisfying ([3]). Also, Cij k + Cj ki + C ki j = 
\{Bij k + Bj ik ) + \{Bj ki + B k ji) + \{B ki j + B ik j) = ^(Bij k + Bj ik - Bj ik + B k ji~ B k ji- Bij k ) = 0, 
satisfying (gj). Hence, ip(B(V)) cC(V). 

We now show C(V) c Choose any C = Cjj&e 1 ® e J ® e fc e C(V) such that Cy* 

satisfies ([3]) and It suffices to choose a tensor i3 = Bij k e l ® e- 7 ® e fc e £>(F) such that 
ip(B) = C. Let Sjjj = for all i. For all i + k, let Ba k = Cu k , B^i = -Ca k and B k a = 0. For 
all i < j < k, let 

Bij k ~ ^Cijk, Bfcij = —ICjki, Bjki = 0, 
Bikj — ~'2>Cijk, Bfcji — 2Cjj : i, Bjik — 0. 

Then, it is easy to show ip(B) = C for the -B constructed above. Hence, C(V) <= 
Since c C(V) and C(V) c ip(B(V)), it follows that r/>(B(V)) = C(F). 

Since ^ is onto, dimker^ = dim.B(V) -dimC(F) = 1 )( n 2 ) Lemmadl □ 

For a manifold M, B{M) and C(M) denote the sets of the tensor fields on M such that 
at each point q e M, B(M) q and C(M) 9 are equal to B(T q M) and C(T q M), respectively, 
where T g M denotes the tangent space to M at g. 

Lemma 3. 1. For a (0,s)-tensor T on a vector space V , T(v, ... ,v) = /or all v € V if 
and only i/Sym(T) = 0. 

Suppose that a (0,3)-tensor C on V satisfies fi]J. Then, Sym(C) - if and only if 
C satisfies fljj. 

Proof, (a) Trivial by definition of Sym. 

(b) Let C = Cij k e % ® e j ® e fc . Then, = Sym(C) = |(Cyfc + C jfci + C feii + C jifc + C fcii + 
Cifci) = \{Cijk + Cjki + C k ij) o C satisfies ©. □ 

Lemma 4. Zef V = V © V", and S = Sijke 1 ® e J ® e fc 6e a (0, 3) -tensor on V that satisfies 

S(u,v,w) - S(v,u,w) Vu,v,weV. (6) 

Then, the following are equivalent: 
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(a) There exists a (0,3)-tensor C = CV^e* <8> e? <8> e fc i/iai satisfies (OP, (0) a^rf 



C(wi,w 2 ,w) = 5 , (wi,w 2 ,'v) V«i,u 2 e V,iJ e V". 



(7) 



(b) 5(ui,u 2 ,u 3 ) + 5(u 2 ,u 3 ,wi) + S(v 3 ,vi,v 2 ) = Vui,u 2 ,U3 e V"- 

Proof, (a) => (b): Trivial. 

(b) => (a): We prove this in coordinates since it will be conveniently used later. Let 
dimV = re, dim^ = n — m and dimU = m. Use the following three groups of indices: 
i,j,k=l,...,n; a, /3,7 = 1, . . . ,n-m; and a,b,c- (n—m+ 1), . . . ,n. Suppose V = spanjej}, 
V = spanjea} and V = span{e a }. Define as follows: 



Cabc = 0. 



(8) 
(9) 

(10) 
(11) 



It is then easy to show that Cijk satisfies ([3]), Q and (|7j). Instead of C a bc = 0, one can 
alternatively choose arbitrary C a bc such that C a bc - Cbac and C a bc + Cb ca + C ca b - 0. □ 

2.2 New Form of IDA-PBC for Mechanical Systems 

Let {ei, . . . , e n } be the standard basis of M. n and {e 1 , . . . , e n } its dual basis. For convenience, 
we identify W 1 with its dual space (M n )*, but we follow the convention that position vector 
q = q % ei is in M n and momentum vector p = p^e 1 is in (M n )*. For the sake of simplicity we 
assume that every function is smooth. 

Let us consider a controlled Hamiltonian system of the form 







In 




d q H 
d p H 







(12) 



where q,p e M n , u e M. m , m < n, I n is the n x n identity matrix, G{q) is an n x m matrix 
with rankG(^) = m for each q, and 

H{q,p) = \p T M-\q)p + V{q) 



is the Hamiltonian function of the system, where M(q) is an nxn positive definite symmet- 
ric matrix and V(q) is a function called the potential function of the system. We say that 
this system has n degrees of freedom and (n - m) degrees of underactuation. Throughout 
this paper, we assume that (q,p) = (0, 0) e M. n xM n be an equilibrium point of the system to 
be stabilized. The controlled Hamiltonian system (112j) shall be called linear if the matrices 
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M (q) and G{q) are constant and the function V{q) is a quadratic function of q. Otherwise, 
it shall be called nonlinear. 

To achieve stabilization of the equilibrium point at the origin, we have the objective of 
finding a feedback law u = u(q,p) that transforms the system f)12|) to the following desired 
form: 



P 



M-\q)M(q) 
-M(q)M-\q) -G{q)K v {q)G 7 \q) 



d q H 
d p H 







M~ l p j M~ l p j C 



where 



H(q,p) = \p T M-\q)p+V(q) 



(13) 



(14) 



is a desired Hamiltonian with an n x n positive definite symmetric matrix M(q) and a 
function V{q) having a non-degenerate minimum at q = 0; K v (q) is an m x m positive 
definite symmetric matrix, and C = Cijk(q)e l ®e J ®e fc is a (0,3)-tensor field that pointwise 
satisfies Q and ©, i.e., C e C(JR ri ). In ((T3J), it is understood that 



M _i p —i M _i p j C = Ci jk M ie M^p iPr e k , 



(15) 



where M* 3 denotes the (i,j)-th entry of the inverse matrix M" 1 of M. Without loss of 
generality we have symmetrized CV,/% with respect to its first two indices (property ([T]) or 
d3])) because M~ l p appears quadratically in (fl~5j) . Notice that Sym(C) = by statement 2 
of Lemma [3j The function V is called the potential energy of the system (|13p . 

If the transformation of (112D to fjl3[) is possible, then the closed-loop system is at least 
Lyapunov stable with H as Lyapunov function since 



dH 
dt 



. dH . dH 

^ dq ^ dp 



= -p T M- l GK v G T M- l p + C(M- l p, M^p, M _1 p) 
= -p T M- l GK v G T M- l p < 0, 

where C(M~ 1 p, M~ 1 p,M~ 1 p) = by ([2]) or statement 1 of Lemma El Since the term 
M~ l p J M~ l p J C does not change the Hamiltonian H, it is called the gyroscopic force. 



2.3 Comparison of the New Form of IDA-PBC with the Conventional 
One for Mechanical Systems 

In the conventional method of IDA-PBC [H [18], for a given system of the form (|12h . the 
following desired form of controlled Hamiltonian system in place of (|13|) is used 



P 



M- l (q)M(q) 
-M{q)M- 1 {q) J(q,p) - G(q)K v (q)G T (q) 



d q H 

d p H A 



(16) 
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with the Hamiltonian (|14p . where J(q,p) is an n x n skew-symmetric matrix that is linear 
in p. In other words, the (i, j)-th entry Jij(q,p) of J(q,p) is written as 

Jij(q,p) = Jij(q)Pk 

where </*;• satisfies 

k k 

due to the skew symmetry assumption on J(q,p). Define B = Bij k {q)e l <8> e? ® e k e i3(M n ) 
by ' 

-Bfcij = JjiMik, (18) 

which is indeed in i3(M n ) by (]17p . Since M is invertible, the relationship in (I18p between 
(Bijk) and (Jy) is bijective. 

Let -F"' = -F/e 2 be the force term J(q,p)d p H in (|16p . which can be written component- 
wise as follows: 

if = Jtffi'pwt 

= B kji M ks W l p s p, 

= i(S fci< + B jki )M ks M je PsPe , (19) 

where the symmetrization in (I19D with respect to indices k and j is valid since Ff is a 
quadratic function of M~ l p. Observe in (|19p that it is the symmetric part ^(B k ji + Bj^) of 
Bkji with respect to the first two indices that essentially contributes to the dynamics (fT6l) . 
but the skew-symmetric part \{B k ji - Bj^) of B^ with respect to the first two indices 
do not appear in the dynamics at all. Since (-Bjjfc) and (Jy) are isomorphically related by 
(fTHj) . it follows that only part of (Jy) contributes to the dynamics (fT6]h 

The natural question that now arises is how to express the essential part of («/£•) 
that appears in the dynamics (|16p . removing all the unnecessary components. Define 
C = C^e* ® e J ' <8> e fc e C(R n ) by 

C(n, w, w) = -(B(u, v, w) + B(v, u, w)) Vu, t> , to e R n , 

or in coordinates 

Cyfc = ~(Bijk + #jifc), (20) 

where 5 = B ijA e* ® e J ' <g> e k e i3(M n ) is the (0, 3)-tensor field defined in (USD. By Lemma 
the tensor field C defined above indeed belongs to C(W l ), and such a C expresses exactly 
the symmetric part ^(Bj,ji + Bj^) that contributes to the dynamics (|16p . Moreover, from 
Lemma [21 we can see that the use of the gyroscopic force M _1 pjM _1 p_iC with C e C(M. n ) 
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in (|13p contains a smaller number of components than (Jh) but still makes the same 
contribution to the dynamics as the term J(q,p)M~ l p does. Concretely speaking, (Jy) or 
equivalently (B^) contains "("^M" -2 ) unnecessary component^] that contribute no terms 
to the dynamics (|16p. Hence, we shall exclusively use the new form (|13p in the rest of the 
paper. 



2.4 Matching Conditions 

We study the synthesis problem: what are the matching conditions for H and C to satisfy 
in order to transform (I12D to f 1 1 3 [) ? 

Let G x (q) denote an (n - m) x n matrix whose rows span the left annihilatoiH of the 
column space of G(q) such that rankG 1 = n-m and G L G = 0. For the sake of simplicity, 
we also denote by G 1 (q) the left annihilator of G(q). Then, by comparing equations (|12p 
and (|13p and collecting terms of equal degrees in p, we obtain the matching conditions: 

= G L (d q V - MM~ x dqV) (21) 
- G L {d q {p T M' l p) - MM' l d q (p T M' l p) + 2M~ 1 p j M^p j C), (22) 

where the first set is called the potential matching conditions and the second the kinetic 
matching conditions. The difference between the two dynamics (|12p and (|13|) is taken care 
of by the feedback law: 



it 



(G T Gy 1 G T (d q H - MM^dqH - GK v G T d p H + M~ x p j M _1 p j C). (23) 



The number of PDEs in the kinetic matching conditions is n ( w+1 ^" m ) if we simply set 
the coefficients of PiPj's to zero in (|22p . However, we can reduce the number of PDEs in 
the kinetic matching conditions by decomposing them into two groups: one without any 
entries of and the other with some entries of C^. For this purpose, let us introduce 
a (2, l)-tensor field A = A l ^(q)ei <8> e 3 - ® e k and a (0,3)-tensor field S = Sij k (q)e l ® e j ® e fc 
that are defined as follows: 



and 



S(u,v,w) = A(Mu,Mv,w) \/u,v,w€R n (25) 



1 Notice that n ( n 1 K" 2 ) > o if and only if the number of degrees of freedom n is greater than 2. 

2 A left annihilator of a subspace [/ of a vector space V is a space defined by {a 6 V* («,tt} = OVtt e [/}. 
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or in coordinates 



Sijk = A r k s M ri M S j 

1 rr ,MdM rs ^~ -~~ ldM 



Notice that 



1— H dMa ldM rs , . 

= - 5 M « M "a-r - 2^ M " M ^'- (26) 

= j4^, and - Sji^. 

Then, the kinetic matching conditions in (|22p can be written in the following compact 
form: 

C(M' l p, M _1 p, w) = A(p,p, w) Vp £l",fflf G 1 
-» G(u, u, w) = Mm, w) Vu et n ,w£ G 1 

^C(u,u,to) =S(u,u,w) Vu6M",w6G i (27) 
o C(«, u, w) = S(u, v, w) Vn, v e M n , 10 e G 1 (28) 

where the equivalence in (|27|) comes from (|25p. and the equivalence in (|28p is due to the 
symmetry in their first two indices of both G and S and the polarization technique, i.e, 
2C(u, v,w) = C{u + v,u + v,w) - C{u - v, u - v, w). The kinetic matching conditions (|28|) 
can be regarded as equations defining the tensor field G on W 1 x 1" x G 1 in terms of M 
and M. By Lemma HI there exists a (0, 3)-tensor field G = Cijk(q)e l <8> e J ® e fc satisfying 
(P), © and (J28D if and only if S satisfies 

S(u, v, w) + S(v, w, u) + S(w, u, v) = Vit, «, w e G 1 , (29) 

which is a set of PDEs for M, not containing any components of G. The number of PDEs 
in ([29D is 

(n - m + 2){n - m + - m) 
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(30) 



which is smaller than w ( w+1 K n m > ; the number of PDEs in (|22p before their re-grouping. 
Once the PDEs (|29p are solved for M, then the tensor field G can be determined purely 
algebraically, following the procedure in the proof of Lemma HI 

This approach is an improvement of the conventional IDA-PBC method in [HdH]. In 
the earlier works, the decomposition of the kinetic energy into the two groups was never 
done, so they often worked with the n ( n+1 )( n m ) PDEs in (|22p regarding J(q,p) in (|16p (or 
equivalently G in the new form (|13p ) as free parameter. For the purpose of comparison, let 
us consider the case of one degree of underactuation, i.e., n — m = 1. In the conventional 
IDA-PBC procedure there are normally more than one PDEs for M coming from the 
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kinetic matching conditions ([22]); see pQ for example. In our method, however, the number 
of PDEs for M is always one by (|3U|) . For example, if G(q) is spanned by {e 2 ,... ,e n }, 
then our kinetic matching condition in (|29p becomes 



MuM*^ + ^M rl M sl = 0. (31) 

This is only one quasi- linear PDE for only M\\. Compare this with the result in [I], where 
the number of kinetic matching PDEs increases as the number of underactuation degree 
increases. This shows the superiority of the new reduced kinetic matching conditions (|29p 
that comes from the decomposition of the kinetic matching condition. 

From the discussions in £ 12.31 it follows that (|29|) could be also obtained with the con- 
ventional form of IDA-PBC (I16j) . but it has never been done anyway. Even with this 
conventional approach, one would have to go through the quadratic gyroscopic force via 
equations (fl~8l) - ([20]) and Lemmas [2] and H] directly or indirectly. By introducing the 
quadratic gyroscopic force from the outset, we avoid this unnecessary detour. 

2.5 Matching Conditions for Integrable G(q) 

In §2.41 we studied how to reduce the number of PDEs in the kinetic matching condition, 
by eliminating the gyroscopic term C. We now show that the number of partial derivatives 
of Mij that appear in the kinetic matching PDEs, can be reduced when the co-distribution 
generated by the column vectors of G(q) in the dynamics (|12[) is integrable, by which we 
mean the integrability of G{q). 

Suppose that G(q) = span{e n ~ m+1 , . . . , e n } for each q such that G L (q) = spanjei, . . . , e n - m }. 
Let us use the following three groups of indices: i,j,k,... = 1, . . . ,n; a, /3,7, . . . = 1, . . . , n-m; 
and a,b,c, . . . = n — m + 1, . . . , n. Then, the potential matching PDEs (|2ip are written as 

dV • • dV 

^--M m M iJ ^- = 0. (32) 
dq a dqJ V ; 

The new kinetic matching PDEs (I29D are written as 

Sap-y + Sfrya + S^aft = (33) 

where 

<? l M M adMa ? ldMrS M M 

Notice that the partial derivatives of only M a p appear in this set of kinetic matching PDEs 
([33]) whereas the partial derivatives of M aa or M a b do not appear. It is also remarkable 
that M a ;,'s do not appear at all in the matching PDEs in (|32|) or (|33|) . 

We now summarize the design procedure for IDA-PBC in the case when G{q) = 
span{e n - m+1 ,...,e n }. 
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PI. Solve the matching PDEs in (J32J) and for M(q) and such that M(q) is 

positive definite and V(q) obtains a non-degenerate minimum at q = 0. (For local 
positive definiteness of M(q), it suffices to have positive definiteness of M(0)) 

P2. Compute C following © - using S ijk in (1261) . 

P3. Choose an arbitrary m x m positive definite symmetric matrix 

P4. Compute the feedback law ([23)) . 

This procedure can be easily adapted to the case of non-integrable G(q). 

2.6 Case of Underact uat ion Degree One 

We consider the following problem: when does the design procedure in §2.51 produce a 
(locally) positive definite symmetric matrix M and a potential function V having a non- 
degenerate minimum at q = 0? In general this is a difficult problem. However, complete 
answers are available for two cases: 1. when the controlled Hamiltonian is linear and 2. 
when the degree of underactuation is one, i.e, n-m- 1. The linear case is important not 
only by itself but also for the nonlinear case because the linear result is used in choosing 
an initial condition for the matching PDEs for the nonlinear case so that H obtains a non- 
degenerate minimum at the equilibrium point. If the degree of underactuation is more than 
one, then solving the matching PDEs becomes a challenging job, which would normally 
involve the operations of prolongation and projection on the system of matching PDEs; 
refer to [19J for definition of prolongation and projection. Luckily, when the degree of 
underactuation is one, no knowledge of formal theory of systems of PDEs such as that in 
|19j is necessary, but a simple application of Frobenius' integrability theorem suffices. 

Let us first take the linear case. Recall that the controlled Hamiltonian system (|12p is 
called linear if M is a constant matrix, V is a quadratic function of q, and G is a constant 
matrix. For IDA-PBC of linear controlled Hamiltonian systems, the gyroscopic term C in 
(|13p is not necessary since it would produce nonlinear terms in the dynamics. 

Lemma 5. Suppose that the controlled Hamiltonian system is linear and has arbitrary 
degrees of underactuation. Then, it can be transformed to the desired form $13\) with a 
constant positive definite symmetric matrix M, a positive definite quadratic function V , 
and C = if and only if the system |Up is controllable or its uncontrollable dynamics are 
oscillatory^. Moreover, the closed-loop system < ti3J) is exponentially stable if and only if the 
original system is controllable. 

Proof. The Lagrangian equivalent of this theorem is proved in [9j, and it can be easily 
adapted to the Hamiltonian case. See also [12] for a brief proof. □ 

3 A linear dynamics x - Ax is called oscillatory if A is diagonalizable and each eigenvalue of A is a 
non-zero purely imaginary number. 
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We now consider the case of nonlinear controlled Hamiltonian systems with underac- 
tuation degree one. 



Theorem 1. Consider the controlled Hamiltonian system hi Sty with one degree of under- 
actuation. Let be its linearization at the equilibrium (q,p) = (0,0). The system can 
be transformed to the desired form \13\) with H having a non- degenerate local minimum 
at (q,p) = (0,0) ; K v = Kj being positive definite, and C e C{W l ), if and only if T, e is 
controllable or the uncontrollable dynamics of T, £ are oscillatory. Moreover, the closed-loop 
system \13\) is (locally) exponentially stable if and only ifTr is controllable. 

Proof. For the sake of simplicity of presentation, we assume that G{q) is integrable such 
that G(q) = { e 2 , . . . , e n }. Then, the potential and kinetic matching conditions in ()32p and 
can be written as 

. . dV dV 

j'- M 'w-w (34) 

+ ^B a M H - 0. (35) 

Write the Hamiltonian H of the linearization of (|12|) at the origin as 

H l = ^p T M- l {Q)p+^q T D 2 V{Q)q 

where M _1 (0) is the evaluation of M~ 1 (q) at q = 0, and D 2 V is the second-derivative 
matrix of V. Suppose that is controllable or its uncontrollable dynamics are oscillatory. 
Then, by Lemma [5] there are two constant positive definite symmetric matrices M and S 
such that 

d 2 v 

dq 1 dq k 



M li M ij (0)S jk -^- r ^(0) = (36) 



which is the potential matching condition for T, e . The kinetic matching condition trivially 
holds for £ since M and M are constant matrices. Notice that the set of PDEs in (plj) and 
(I35p for V and M\\ is integrable, which can be easily checked by applying Frobenius' the- 
orem; refer to |17] for Frobenius' theorem. We now impose the following initial conditions 
on M and at q = 0: 

M(0) = M, dP(0) = 0, D 2 V(0) = S. (37) 

The initial conditions at q = are compatible with the PDE (|34p since the differentiation 
of (|34p and its evaluation at q = yields 

OqWq K oq^oq^ 
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where dV(0) - is used. Equation ([38]) has the same structure as ([36]) . Hence, by 
Frobenius' theorem there exist M and V that satisfy (|34p . (|35j) and (|37|) . Due to the initial 
conditions (|37p . the desired Hamiltonian = ^p T M^p + V(q) obtains a non-degenerate 
local minimum value at (q,p) = (0,0). Choose C, K v and u following the procedure in 
We have successfully transformed the given system (|12p to the desired form ()13[) . 
Suppose now that the given nonlinear controlled Hamiltonian system fjl2[) can be trans- 



formed to the desired form (|13p with H = \p T M 1 p + V(q). Through linearization, we can 



2 

see that the linearization S of (|12p gets transformed to the following linear dynamics: 



P 



M- 1 (0)M(0) 
-M(0)Ar 1 (0) -G(0)i^(0)G T (0) 



(39) 



where H l = \p T M~ 1 ({))p+ ±q T D 2 V(0)q with both M(0) and D 2 V(0) positive definite 
and symmetric. Hence, the linear controlled Hamiltonian system is controllable or its 
uncontrollable dynamics are oscillatory by Lemma [5j This completes the proof of the first 
statement that the system (I12D can be transformed to the desired form (I13p if and only if 

is controllable or the uncontrollable dynamics of are oscillatory. 

We now prove the statement on exponential stability. By the Lyapunov linearization 
theorem, the closed-loop system (|13p is exponentially stable if and only if its linearized 
dynamics, say the one in (I39p . is exponentially stable. From Lemma [5j we know that (I39p 
is exponentially stable if and only if the linearization of the original nonlinear controlled 
Hamiltonian system (112p is controllable. This completes the proof. 

□ 

Theorem [1] is important because it gives necessary and sufficient conditions for Lya- 
punov/exponential stabilizability by IDA-PBC for the case of underactuation degree one. 
Moreover, the conditions therein are on the given system (|12() that can be easily verified in 
advance before any attempt to find a desired system (|13[) . In the past, they had only suffi- 
cient conditions for stabilizability by IDA-PBC, so they were inconclusive on stabilizability 
by IDA-PBC when the sufficient conditions did not hold. 

Theorem [T] applies to a wide range of systems including the inverted pendulum on a car, 
the Pendubot, the Furuta pendulum, the ball and beam system, and the planar vertical 
take off and landing aircraft, the linearization of each of which is controllable. We remark 
that the resultant feedback controller (I23p is nonlinear. 

2.7 Examples 

We consider the inverted pendulum on a cart. The mass matrix M and the potential 
function V are given by 

M _\ M x l 2 M^cosq 1 ' 
Milcosq 1 Mi + M 2 
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and 



V{q) = M 1 gl{l-cosq 1 ). 
Since control is given in the q 2 direction, we have 



G(q) = 



Our main purpose is to illustrate the design procedure in §2.51 s ° we choose the following 
values of parameters: 

Mi = M 2 =£ = l, 5 = 10. 

One can easily check that the linearization of this system at the origin is controllable. 
Hence, by Theorem [T] we can exponentially stabilize this system with the method of IDA- 
PBC. 

Let us now construct a control law following the procedure in §2.51 The potential and 
kinetic matching conditions in (j32H and (I33p are given, after simplification, as 

Qy Qy 

(2M n - Mm cosg 1 )— + (-Mn cosg 1 + M 12 )— = 10 sing 1 (2 - cos 2 q 1 ), 
oq l oq l 



and 



(2M n - Mi 2 cos q 1 )^- + (-M n cosg 1 + M 12 )^i 
oq L oq L 

_ 2smq 1 (2M n - M 12 cos q 1 )(M 12 -M u cos q 1 ) 

2 - cos 2 q 1 



A solution is found to be 



Mn = 2cosV -e, Ml 2 = (4-e)cosg\ 

2 



V 



10 



cos q + I q 



2 sing 1 



where < e < 2 is a constant. To make M positive definite at least locally around q = 0, we 
choose an M22 such that M 22 > (Mi2) 2 /Mn at q = 0. For example, 



(4-e) 2 cos 2 q 1 
2 cos 2 q 1 - e 



M 22 = K + 



where K > is a constant parameter. It is easy to see that (q,p) = (0,0) is the minimum 



point of the Hamiltonian H(q,p) = -^p M p+ V(q) over the set /, 



£ 2 , where 



L ■- cos 
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Notice that I e -> (-■ | , ^) as e -> 0. Compute the gyroscopic tensor C using ([8]) - (fTTj) and 
(I26p as follows: 



Cm 


= C222 


Cm 


= C211 


C221 


= 5*221 


C112 


= -25] 


C122 


= C212 



s 



121 



121 
= -^221, 



where 



5*121 
5*221 



M lr M T 



-~\MvrM rs 



dM 12 


dM rs 


dq s 


+ Oq 1 


dM 22 


dM rs 


dq s 


+ dq 1 



— M ri M s2 J 



—M r2 M s2 



Choose a If t, = i^J > and compute the feedback control law (|23|) . The closed- loop system 
is exponentially stable by Theorem [TJ 

Let us consider another example. Suppose that we are given a controlled Hamiltonian 
system of the form (|12[) with 



M(q) 



5 + cost/ 3 sin(g 1 -g 2 ) sm(q 3 - q 1 ) 
sin(g 1 -g 2 ) 5 + cos(g 1 -g 3 ) sing 2 



sin(g 3 - q 1 ) 



smg 



5 + cos(g ) 



V(q) = cos(g 1 + q 2 ) + cos(g 2 + q 3 ) + cos g 3 , 



and 



G(q) 



smq 



smq 



sin g 
1 



where g = (q l ,q 2 ,q 3 ) e M 3 . We want to stabilize the equilibrium at the origin in IR 3 x R 3 . 
This system has one degree of underactuation, and its linearized system at the origin is 
controllable, which is easy to verify. Hence, by Theorem Q] the origin can be exponentially 
stabilized with the IDA-PBC method. 



3 Conclusions 

We have improved the method of IDA-PBC in several ways. First, we have showed that 
there is redundancy in the skew-symmetric interconnection matrix in the conventional form 
of IDA-PBC, and then have replaced the skew-symmetric interconnection matrix term with 
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a gyroscopic force to remove the redundancy. We have derived the matching conditions 
and decomposed them into two parts. As a result of the decomposition, we have come up 
with a smaller number of kinetic matching PDEs than those in the literature. Moreover, 
a smaller number of entries of the desired mass matrix appear in the kinetic matching 
PDEs when controls are given in coordinate directions. Easily verifiable necessary and 
sufficient conditions are given for Lyapunov/exponential stabilizability by IDA-PBC for 
all linear controlled Hamiltonian (or simple mechanical) systems with arbitrary degree of 
underactuation and for all nonlinear controlled Hamiltonian (or simple mechanical) systems 
with one degree of underactuation. A synthesis procedure with IDA-PBC was provided 
and illustrated on the system of an inverted pendulum on a cart. We plan to study the 
case of higher degree of underactuation and force shaping for IDA-PBC in the future. 
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